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We investigate the entanglement entropy (EE) of gapped S = 1 and S = 1/2 spin chains with 
dimerization. We find that the effective boundary degrees of freedom as edge states contribute 
significantly to the EE. For the S = 1/2 dimerized Heisenberg chain, the EE of the sufficiently 
long chain is essentially explained by the localized S = 1/2 effective spins on the boundaries. 
As for S = 1, the effective spins are also S = 1/2 causing a Kennedy triplet that yields a lower 
bound for the EE. In this case, the residual entanglement reduces substantially by a continuous 
deformation of the Heisenberg model to that of the AKLT Hamiltonian. 
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1. Introduction 

The concept of entanglement in quantum mechanics 
has significant importance in quantum information tech- 
nologies, where quantum resources are used for process- 
ing information in novel ways. Recently, it has been clar- 
ified that the entanglement of a quantum state is one of 
the most important properties not only in quantum in- 
formation science but also in condensed matter physics. 
Entanglement entropy (EE) measures the entanglement 
of a single quantum state, which is defined as the von 
Neumann entropy of a reduced density matrix of a sub- 
system. EE has been applied to a condensed matter sys- 
tem in order to understand quantum phases. There are 
several works 1 4 that have reported the EE near or at 
a quantum critical point. As for the system size depen- 
dence, the EE of gapped spin chains saturates, although 
that of critical spin chains shows a logarithmic diver- 
gence. Recently, the EE was used to detect the quan- 
tum dimension, which is a property of the topological 
order. 5 7 The topological order is a new order that can- 
not be characterized by a classical local order parameter 
and a symmetry breaking by the conventional Ginzburg- 
Landau theory. 8 Recently, one of the authors proposed 
the use of quantum objects such as quantized Berry 
phases to define the local order parameters that do not 
require classical symmetry breaking since they are gauge- 
dependent quantities unlike the classical ones. 9 12 

It is also interesting to study the EE from the view- 
point of bulk-edge correspondence. Although the EE is 
a physical quantity of the bulk, there is a relation be- 
tween the EE and the edge states that appear in a 
topologically nontrivial system with boundaries. There 
are several works on the bulk-edge correspondence in 
many kinds of physical systems such as the 2D quan- 
tum Hall effect, 13,14 polyacetylene, 15, 16 and the Haldane 
spin chain. 17 ' 18 The EE and the Berry phase of a topo- 
logically nontrivial system are closely related through the 
edge states. 7 
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In this paper, we investigate the relation between the 
saturated EE in a gapped system and the corresponding 
edge state through the bulk-edge correspondence. The 
systems chosen are three different types of gapped quan- 
tum spin chains: (i) the 5=1/2 dimerized Heisenberg 
chain, (ii) the S = 1 XXZ Hamiltonian with the on- 
site anisotropic term D J^Li (^f) 2 \ an d (hi) the Affleck 
Kennedy Lieb Tasaki (AKLT) model. 19 

The S = 1 XXZ chain near the Heisenberg point is a 
typical Haldane spin chain, which has a finite gap in an 
appropriate range of the parameters. There are several 
other phases such as the "large- D phase" and the "Neel 
phase." They are characterized by the hidden Z2 x Z2 
symmetry breaking, which is nonlocal. 20 ' 21 The string 
order parameter is a nonlocal order parameter; 22 there- 
fore, it is not a conventional classical order parameter. 
We see that the EE reflects these kinds of hidden symme- 
try breaking, which relates to the structure of the edge 
states. 

2. EE of gapped system and edge state 

The EE (the von- Neumann entropy) used in this paper 
is defined as follows: 

S = - (log p A ) = -Tr A [p A log p A ], (1) 

p A = Tvb[p) , where p is a density matrix of the pure state 
We assume that the system consists of subsystems 
A and B. The reduced density matrix pA is obtained 
by tracing out the degrees of freedom in subsystem B 
from the total density matrix p. The EE quantifies the 
information about how much the state is entangled be- 
tween the subsystems A and B. We consider to be 
the ground state (GS), assuming that it is unique. 

We divide a one-dimensional chain into subsystems A 
and B (Fig. 2) that denote the upper and lower parts, 
respectively. Since taking a partial trace over a spatially 
separated subsystem results in effective boundaries, it 
induces effective localized degrees of freedom as edge 
states. Then, we evaluate the EE as S ~ Sbulk + S e d ge 
where Sbulk an d S e d ge denote the contribution from the 
bulk and edge, respectively. This is because the degrees of 
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Fig. 1. The EE of the periodic dimerized Heisenberg chains. J\ = 
sin and J2 = cos 0, and the system size is TV = 8, 12, 16, 20. The 
base of the logarithm is 2. 
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Fig. 2. (Color Online) Schematic diagram of taking a partial 
trace in the periodic dimerized Heisenberg system. The degrees of 
freedom are revealed on the edges on which the effective S = 1/2 
spins localize. 



freedom originate not only from the edges but also from 
the bulk part of the system. We hypothesize Sbulk > 
and, we obtain the relation S > S e d ge ^ which means that 
the EE has lower bound of S e d ge - The contribution from 
the edges becomes important in a gapped system, while 
that from the bulk is essential in a gapless system be- 
cause of the infinite correlation length. 7, 23 

3. Evaluation of EE for several spin chains 

We calculate the EE for the three spin chains and show 
that the edge state picture of the EE is valid for these 
systems. We obtain the ground state of the systems nu- 
merically using the Lanczos method in an invariant sub- 
space with S* otal = Sf = 0, where Sf denotes the 
z-component of the spin operator on the i-th site. The 
EE is evaluated by diagonalizing pA using the House- 
holder method. We calculate the logarithm in base 2. 
We impose periodic boundary conditions for the systems. 
These models exhibit various behaviors that can be dis- 
tinguished as different quantum phases by the EE. 

5 = 1/2 Dimerized Heisenberg Model: The 
Hamiltonian of the 5 
system is given by 



1/2 dimerized Heisenberg 



H 



S=1/2,D 



N/2 

£ 

2=1 



(J\S 2 i-l ■ S 2 i + J2S21 ' «2i+l) ,( 2 ) 



where Si are the 5=1/2 spin operators. J\ and J2 are 
parametrized as J\ = sin# and J2 = cos#, respectively, 
where G 5 1 . We consider the case of — tt/2 < 9 < tt/2 
in this paper. The ground state is unique in any 6. 26 In 



the case of < 8 <C \ , the ground state is composed of 
an ensemble of N/2 singlet pairs and the energy gap is 
finite (~ The system is equivalent to the isotropic 

antiferromagnetic Heisenberg chain at 6 = j. This is the 
only gapless point. In the ferromagnetical strong cou- 
pling, 6 — > — 77 + 0, the system is effectively considered 
as the 5 = 1 Heisenberg chains. 26 

We now discuss the EE of the system. Fig.l shows 
the 6 dependence of the EE. Fig. 2 shows the subsystems 
considered in this study. As shown in Fig. 2, a partial 
trace is performed across two J 2 coupling bonds. The 
EE tends to diverge as the system size increases at the 
critical point 6 = \ (isotropic Heisenberg point). This 
behavior of the EE has been discovered by Vidal et al. 3 
One can clearly observe the EE of the edge states, 2 log 2, 
near ~ 0, while the EE converges to in the region 
of — > f . In the former limit, the system becomes a 
collection of N/2 singlet pairs with two almost free spins 
near the two boundaries. Therefore, the EE in this region 
is given by that of the two free spins near the boundaries 
2 log 2. In the latter limit, the system also becomes a 
collection of N/2 singlets; however, there are no free spins 
in this case. Then, the EE vanishes. 

These results can be interpreted as follows. The EE 
reflects the local degrees of freedom that appear near 
the boundaries when the system gets truncated. Then, 
we speculate that a lower bound is obtained for EE as 

5 > logg, where g denotes the above mensioned degrees 
of freedom near the edges. In this case, this local degree 
of freedom arises due to the two 5 = 1/2 spins localized 
on the edges as g = 4 as \q)l ® \P)r (a, P =T> I)- 

The low energy behavior of this system is identical to 
that of the 5 = 1 Heisenberg model in the region of 

6 — > — I + 0. 26 We expect that the edge states of the 
5 = 1 Haldane phase are similar to those of the 5 = 1/2 
dimer phase as the EE is larger than 2 log 2, which gives 
the lower bound. 

S=l XXZ Chain with D-term: The Hamiltonian 
of the 5 = 1 XXZ chains with an on-site anisotropic term 
is given by 



N 

H S =1 = J ^ { S i S i+l 
1 



qy qy 



A5f 5f +1 



N 



(3) 



where Si are the 5=1 spin operators. The energy gap as 
a Haldane gap is finite near the point A = 1 and D = 0. 2T 
There are several phases known as the "Haldane phase," 
the large-D phase and the Neel phase in the range of 
A > with a gapped-gapped transition (Gaussian tran- 
sition) between the Haldane phase and the large-D phase 
and the Ising transition between the Neel phase and the 
Haldane phase. These three phases are distinguished by 
local and nonlocal order parameters 22 or by the breaking 
of the hidden Z2 x Z2 symmetry that appears on per- 
forming a non-local unitary transformation. 20 ' 21 Quan- 
tum order parameters to characterize the states are also 
defined. 9 11 A system with an open boundary condition 
has the localized effective spin- 1/2 edge state in the Hal- 
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Fig. 3. (Color Online) (a) The EE of the S = 1 XXZ chain with 
the D term (N=10). A represents the anisotropies along the S z 
direction. D is the coefficient of the D term. The base of the 
logarithm is 2. (b) Cross section of Fig. 3 along line A. A = 1.0, 
N=4,6,8,10,12.(c) Cross section of Fig.3 along line B. D = 0.0, 
N=4,6,8,10,12. 



dane phase, 17, 18 and the low energy effective Hamilto- 
nian can be written by an effective coupling of the two 
S = 1/2 edge spins. 17 

Fig. 3 (a) shows the EE of the XXZ spin chain with the 
D term, which is calculated in the subspace Sl° tal = 
by changing the anisotropies A and the coefficient of the 
D term D with a system size N = 10. Figs. 3(b) and (c) 
show the cross section of Fig.3 (a) along the correspond- 
ing planes (indicated by lines A and B in Fig. 3 (a)) and 
the system size dependence of the EE. 

From Fig. 3 (a), it can be observed that there are 
three regions. The EEs of the large-D phase, the Neel 
phase, and Haldane phase are approximately 0, log 2, 
and 2 log 2, respectively. We now discuss the relation be- 



tween the EE and each phase. We use bases S£\±)i = 
±\±)i, Sf\0)i = 0. The large-D phase (blue region) 
has a "quantum-gas-like" GS. In the large- D limit, the 
quantum fluctuation ^ [S^S^ + S^ S^] can be ig- 
nored. Then, this state can be approximately described 
by |large-D) = ®f =1 |0)i, which results in a vanishing 
EE because this state is the direct product state. This 
state does not have any degrees of freedom in the sub- 
system created by taking a partial trace. In the Neel 
phase, the GS is Neel ordered as |Neel) = |(-l)*)i 
(quantum-solid-like). We consider the degrees of freedom 
of the edge states created by taking a partial trace. We 
denote the spin on the left and right edge as Sl and Sr, 
respectively. Then, Sl is correlated with Sr in this state 
due to the Neel order. This implies that there are two 
degrees of freedom. Then, the EE is described by the ef- 
fective degrees of freedom as ® |— )r, ® 
(subsystem size : even) and \—) L (g) \-)r, ® 
(subsystem size : odd). This yields logg (g = 2) as the 
EE. Finally, the Haldane phase is a quantum liquid. The 
states which include +, — particles as many as the num- 
ber of mainly contribute to GS. There are localized 
effective spin- 1/2 edge states in the Haldane phase. The 
EE is larger than 2 log 2, which can be interpreted as the 
contribution from the edge state. The rest of the contri- 
bution 5 — 2 log 2 comes from the bulk. This interpre- 
tation is consistent with the result. The EE reflects the 
spontaneous breaking of the hidden Z2 x Z2 symmetry, 
which plays an essential role in Haldane spin chains 20 ' 21 
and it also relates to the edge states. In summary, the 
EE has a lower bound of logg, where g is the number 
of elements in the group of spontaneously broken sym- 
metry (g = 1 for the large-D phase, g = 2 for the Neel 
phase, and g = 2 2 for the Haldane phase), at least in the 
calculated regions. 

We notice that the EE is lower than 2 log 2 in the Hal- 
dane phase when total system size is small, as shown 
in figs. 3 (b) and (c). This low value of the EE is due 
to the correlation length of the system and the system 
size. When the correlation length is long as compared 
to the total system size, Sl and Sr are effectively cou- 
pled. Then, the degrees of freedom on the edges decrease 
because both the sides of the edge states are not inde- 
pendent in the presence of substantial coupling. 

As discussed above, the bulk-edge correspondence is 
confirmed through the trace-out operation, and the EE 
shows the degrees of freedom in the subsystem. An edge 
state in the open boundary condition is one of the impor- 
tant characteristics of the Haldane materials. The above- 
mentioned phases of the periodic spin chains are well 
characterized by the EE that reflects the edge state of 
each corresponding open spin chain. We found that the 
bulk and edge are inextricably linked in these kinds of 
spin chains. We confirmed that the Haldane phase can be 
distinguished by the EE intermediating the information 
of the edge states. 
The S = 1 Heisenberg Model to the AKLT Model: 
We also calculate the EE of the AKLT model numerically 
by modifying the parameter (3 (Fig. 4). The Hamiltonian 
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is given by 

N 

Haklt = jJ2{ Si ' + /? s * ■ ^+i) 2 ) > ( 4 ) 

2=1 

where are the 5 = 1 spin operators. There are solvable 
points at f3 = l/3 19 and /3 = ±1. 28 32 The ground state 
can be explicitly written as the valence-bond solid (VBS) 
state at the point /3 = 1/3. 

Fig. 4 shows the EE of the AKLT model obtained by 
changing the parameter (3. The result shows that the 
value of EE is minimum at the point (3 = 1/3, at which 
the ground state is the VBS state. The EE is lower than 
2 log 2 when the system size is N = 4, 6, 8, 10. This is due 
to the correlation between the edge spins (Sl and Sr), 
which is also discussed in the previous model. According 
to the edge state picture of the EE, S > 2 log 2 when the 
substantial coupling of the spins vanishes (in the large N 
limit), as speculated from the numerical results; however, 
one cannot observe this relation from the numerical cal- 
culations. To confirm this relation, we evaluate the exact 
form of the EE of the periodic VBS state analytically. 
The EE of the periodic VBS state is obtained analyti- 
cally by using a transfer matrix technique 33 as follows: 



?N,L 



Xa 



-3X A log X A 
l(l-p L )(l 



X B log A B , 

N-L\ 



P 



Ar = 



4 

1(1 



1-p 
3p L )(l 



N-l 



3p 



N-L\ 



1 



P 



N-l 



(5) 
(6) 

(7) 



p = — - , where N and L denote the total system size and 
the length of the subsystem, respectively. It is evident 
that this expression is invariant under the replacement 
U = N — L, i.e., Sn,l = $n,n-l, which is a general 
property of the EE. 

In the limit N — » oo, the EE reduces to the expression 



lim S NjL 

N^oo 



21og2 



-p L )\og(l-p L ) 



1 



(l + 3p L )log(l + 3p L ), 



(8) 



which is exactly the same as the EE of the VBS state 
with two spin- 1 /2s on the boundary. 34 

Next, we consider the thermodynamic limit of Sn,l 
with a = L/N fixed. We obtain limjv^oo Sn,Ncx = 2 log 2 
(a = L/N fixed). We can see that the EE does not di- 
verge. This result is consistent with the equation S > 
2 log 2 as a consequence of the edge state picture. Ac- 
cording to the edge state picture, the EE in the AKLT 
model receives no contribution from the bulk, while the 
contribution from the two edges is present. Therefore, 
this is the point at which the EE is lowest in the Hal- 
dane phase, as long as the edge state picture of the EE 
is valid. 

4. Conclusion 

We have explored new aspects of the bulk-edge corre- 
spondence in spin systems on the basis of quantity ob- 
tained from the quantum information, namely, the entan- 
glement entropy (EE). The tracing out can be interpreted 




Fig. 4. (a) The EE of the periodic S = 1 AKLT model to the 
Heisenberg Chain, (b) Size dependence of the EE at the AKLT 
point. The base of the logarithm is 2. 



as the truncation of the system, and the EE includes the 
information regarding how many degrees of freedom the 
truncated state has. In the case of the S = 1 XXZ spin 
chain with the D term and the AKLT model, it is the 
degrees of freedom of the effective 5=1/2 spins local- 
ized at the edges of the chain, which are known as the 
characteristics of the ground state in the Haldane phase 
with an open boundary condition. Futher, we speculate 
that the EE has a lower bound of logg, reflecting the 
effective boundary degrees of freedom g. 
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